Nano-wires with surface disorder: 
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We investigate electron quantum transport through nano-wires with one-sided surface roughness 
in the presence of a perpendicular magnetic field. Exponentially diverging localization lengths are 
found in the quantum-to-classical crossover regime, controlled by tunneling between regular and 
chaotic regions of the underlying mixed classical phase space. We show that each regular mode 
possesses a well-defined mode-specific localization length. We present analytic estimates of these 
mode localization lengths which agree well with the numerical data. The coupling between regular 
and chaotic regions can be determined by varying the length of the wire leading to intricate structures 
in the transmission probabilities. We explain these structures quantitatively by dynamical tunneling 
in the presence of directed chaos. 
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I. INTRODUCTION 

Deterministic dynamical systems with a mixed regular- 
chaotic phase space give rise to many interesting fea- 
tures that are present in neither the regular nor the 
chaotic limit alone. These special features include, in 
the purely classical limit, the emergence of directed chaos, 
where the chaotic part of phase space exhibits a diffusive 
motion with a drift in a specific direction. Interest- 
ingly enough, such a biased diffusion may be realized in 
a mixed phase space even without external driving, re- 
quiring only broken time- reversal symmetry.'' Quantized 
versions of mixed systems also have interesting proper- 
ties, since different regular and chaotic regions in phase 
space, which are completely separated classically, become 
connected by dynamical tunneling.^ Dynamical tunneling 
can be viewed as a generalization of conventional tun- 
neling through potential barriers to multi-dimensional 
non-separable dynamical systems where the tunneling 
path along a specific "reaction coordinate" is, in gen- 
eral, not well defined. Dynamical tunneling in phase 
space plays an important role in spectral properties and 
transport. ^^^^ Dynamical tunneling rates between regu- 
lar regions are substantially enhanced by the presence of 
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FIG. 1: Wire with one-sided surface disorder in a magnetic 
field B, applied perpendicular to the scattering area. The 
solid line shows a regular skipping trajectory and the dashed 
line an irregular trajectory scattering at the disordered sur- 
face. 
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FIG. 2: (Color online) Poincare sections for (a) right-moving 
{px > 0) and (b) left-moving {px < 0) classical electrons in 
the wire shown in Fig. 1. In (a) a single regular island with 
invariant tori (green) can be identified, which is separated 
from the chaotic region (blue dots) by the outermost torus 
(dashed). In (b), no such island appears. The gray-shaded 
part indicates the y-range affected by disorder, cf. Fig. 3. 



chaotic motion, giving rise to the notion of chaos-assisted 
tunneling. Such dynamical tunneling processes were re- 
cently observed with cold atoms in periodically modu- 
lated optical lattices. 

In the present paper, we study transport in a long 
quantum wire with one-sided surface disorder in the pres- 
ence of a perpendicular homogeneous magnetic field (see 
Fig. 1). This system features coexistent regions of reg- 
ular motion (skipping trajectories) and irregular motion 
induced by scattering at the disordered surface. In the 
corresponding mixed phase space (Fig. 2) we find di- 
rected regular and irregular motion, which are quantum 
mechanically coupled by dynamical tunneling. As a re- 
sult, this scattering system gives rise to many interest- 
ing properties which are clearly reflected in the trans- 
mission through the wire — a quantity which is readily 
accessible in a measurement. Surface disordered wires 
and waveguides have recently received much attention, 
both theoretically^""^^ and experimentally.^^"^® This is 
mainly due to the fact that nano-wires, albeit being con- 
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ceptionally simple, are very rich in their physical proper- 
ties. Nano- wires are now also being realized as graphene 
nano-ribbons, for which surface disorder seems to play 
an even more important role than for conventional (semi- 
conductor based) wires. 

In a previous paper''^ we showed both numerically and 
analytically that by increasing the number of open chan- 
nels N in the wire, or cquivalcntly, by increasing the 
wavenumber k-p, the localization length ^ induced by 
surface disorder increases exponentially. The localiza- 
tion length was related to tunneling of the lowest trans- 
verse wire mode (m = 1) from the regular island to the 
chaotic region in phase space (see Fig. 2). The dramatic 
increase of ^ then follows directly from the exponen- 
tial suppression of the tunneling rates in the semiclas- 
sical limit of large fcp (or small de Broglie wavelength 
Ad). In the present paper we explore the behavior of 
higher transverse modes, m > 1. We give for each of 
these modes m a remarkably accurate analytical esti- 
mate for its specific localization length, ^,„. Further- 
more, all modes in the regular island are effectively cou- 
pled to one another by dynamical tunneling in the pres- 
ence of irregular, yet directed, motion: key is here the 
interplay between directed regular motion (to the right 
in Fig. 1) and counter-moving irregular motion directed 
to the left which gives rise to characteristic structures 
in the mode-specific transmission probabilities Tm of the 
current-transporting regular modes m. These intricate 
structures can be accounted for by a simple scattering 
model incorporating opposing directed regular and irreg- 
ular motions and their coupling by tunneling. 

This paper is organized as follows: In section II we 
briefly review the characteristic classical and quantum 
features of our nano- wire, as induced by its mixed phase 
space. In section III we provide analytic estimates for 
the individual mode localization lengths S^m- Techni- 
cal details of the underlying calculations are deferred 
to appendix A. In section IV we analyze the individual 
modc-to-mode specific transmission probabilities T^.m', 
for which a conceptually simple transport model based 
on the coupling by dynamical tunneling between regular 
and irregular directed flow in opposite directions. The 
paper is rounded off by a summary in section V. 

We use atomic units, but include the constants h — 
iTieS = e = Atteo ~ 1 where instructive. 



II. WIRE WITH SURFACE DISORDER 

A. Classical dynamics 

We consider a 2D wire with one-sided surface disorder 
to which two leads of width W are attached (Fig. 1). 
A homogeneous magnetic field B perpendicular to the 
wire is present throughout the system. We choose the 
magnetic field to be directed in negative z-direction. The 
Hamiltonian H ^ ^(p + -A-)^ -I- V{x, y) is then given in 



Landau gauge (A = By x) by 

H = ]^{p, + Byf + ?j-^V{x,y). (1) 

In the leads (x < and x > LI) 

V{x,y) = Vo[Q{-y) + e{y-W)] , (2a) 

and inside the quantum wire Q < x < LI 

V{x, y) - Vo M-y) + Q{y-W + i^x))] . (2b) 

Vb is taken to be arbitrarily large (Vq oo) to represent 
hard wall boundary conditions. To emulate stochastic 
classical scattering at the upper wire surface we choose 
rjix) to be a random variable that is piecewise constant 
for a fixed interval length I (see Fig. 1). We choose the 
value of 'q[x) to be uniformly randomly distributed in the 
interval 

-5/2<r]{x) <5/2. (3) 

Thus, the wire is assembled from L rectangular elements, 
referred to in the following as modules, with equal width 
Z, but random heights /i, uniformly distributed in the 
interval \W — 5/2, W + 5/2]. In the numerical computa- 
tions we use I = W/5 and 5 = {2/3)W. The Hamiltonian 
Eq. 1 is non-separable for a given realization of disor- 
der. It is therefore quite likely that generally in such 
a system, mixed regular and chaotic motion may ensue. 
In the following, we refer to the irregular motion in the 
wire as "chaotic" , although we will only make use of the 
weaker condition of stochasticity and ergodic coverage of 
the corresponding phase space region due to the random 
disorder. The ensemble-averaged value of the Lyapunov 
exponent does not enter the subsequent analysis. 

The classical motion in the wire proceeds on circular 
arcs characterized by the cyclotron radius Vc — Pf/B 
and guiding center coordinate yo ~ ~Px/B, interrupted 
by elastic refiections on the hard wall boundary. For y^ 
sufficiently small or negative such that yo + '^c < W — 5/2, 
the electron performs regular skipping motion, for which 
2/0 is conserved. These skipping trajectories generate a 
directed, ballistic motion to the right (for the direction of 
the B-field and boundary conditions depicted in Fig. 1). 
For trajectories hitting the upper disordered boundary 
(2/0 + rc> W — 5/2), chaotic motion will develop with an 
average drift to the left, as discussed below. 

To visualize the classical dynamics for electrons with 
Fermi momentum pp we choose as Poincarc section a ver- 
tical cut at the entrance of the wire {x = 0) with periodic 
boundary conditions in the x-direction, ' see Fig. 2. The 
resulting section (y, py) for p^ > shows a large regular 
region with invariant tori corresponding to the skipping 
motion along the lower straight boundary of the wire. 
The section for p,j. < shows only an irregular region, 
as all trajectories with p^ < interact with the upper 
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FIG. 3: (Color online) (a) Nano-wire with the regular trans- 
verse modes Xm{y) m — 4,3,2,1 (green) for kpW/TT = 14.6. 
The gray shaded part indicates the j/-range aflFected by disor- 
der, (b) Poincare-Husimi functions of these modes and their 
quantizing tori. 



boundary. The area A in the Poincare section enclosed 
by a torus is given by 



PFrc 



arccos(l - u) - {1 - v)^\ - (1 - vf , (4) 



where vr^ is the y-position at the top of the correspond- 
ing cyclotron orbit. The size Aiog of the regular island is 
obtained for I'max = (W^ — '5/2)/rc. Outside of the regular 
region the motion appears uniformly chaotic. Hierarchi- 
cal structures of island chains are absent. 



The system displays directed chaos,'' 



i.e., the time 



averaged velocity of almost all classical trajectories in 
the chaotic region of phase space approaches a non-zero 
constant Vch for long times. This chaotic drift motion 
arises here as trajectories in the regular island have a 
non-zero average speed Wreg > directed to the right, 
while the velocity average over the whole phase space 



is exactly zero. A,-, 



^rcg 



^chi'ch = 0. Therefore the 



average velocity of the chaotic part must be directed to 
the left, ■^ch < 0. 



B. Quantum description 

Quantum mechanically, the scattering through the 
wire is described in terms of the modes in the entrance 
and exit lead, respectively, see Fig. 3(a). The low- 
est transverse modes of the incoming scattering wave- 
functions reside inside the regular island (Fig. 3(b)). 
Only their exponential tunneling tail in the harmonic- 
oscillator-like potential (see Eq. 1) touches the upper dis- 
ordered surface at y > W — S/2. These regular modes 
can be semiclassically quantized as'^^'"^'^ 
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1/4) with TO =1,2, 



(5) 



(5) 








\ 




(b) 














« 








m = 4 


1 ** 


2 


\ ■ 





10'^ 



10" 



10« 



L 



FIG. 4: (Color online) (a) Averaged conductance (g) vs length 
L of the wire for rc = and kpW/TT = 14.6. The step-wise 
decrease is accompanied by the disappearance of the regu- 
lar modes (lower set of Husimi distributions) and the flood- 
ing of the island region by the chaotic states (upper set). 
The Poincare-Husimi representations to the left (right) of the 
curve correspond to scattering from left to right (right to 
right) lead, respectively, (b) Transmission 7^ = exp({lnTm)) 
of the incoming mode m vs L. The light vertical lines corre- 
spond to the predicted localization lengths ^,n, Eq. 12. 



where A is the area in the Poincare section enclosed by 
a given torus, see Eq. 4, and A ~ rcA/pp is the area in 
position space enclosed by a segment of a skipping orbit. 

The quantum states of the system can be represented 
on the Poincare section by projecting them onto coher- 
ent states coh(j, (y'), which are localized in the point 



Hm{y,Py} = 



w 



Xrniy') co\y^p^){y') dy' 



(6) 



to obtain a Poincare-Husimi representation. Unlike pre- 
vious studies (see [34,35] and references therein) we use 
the states themselves and not their normal derivative 
function as the section does not run along a border of 
the billiard. The Poincare-Husimi representation for the 
lowest lead modes (Fig. 3(b)) clearly displays the local- 
ization on the quantized tori. 

As this system is a non-separable two-degree of free- 
dom system, the one-dimensional states in the transverse 
direction y do not remain stationary in the diamagnetic 
harmonic-oscillator potential of the Hamiltonian in Eq. 1 
(truncated at y = 0,rj). Due to their extension to the 
rough upper wall, the transverse states behave as an 
equivalent time-dependent one-degree of freedom system 
in the y-coordinate with a time-dependent potential due 
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to stochastic variations of the upper wall. Accordingly, 
each mode has a characteristic decay time 7^ which is 
controlled by the overlap of the tunneling tail with the 
rough wall. Modes with small m live near the center of 
the island and therefore only couple weakly having small 
tunneling rates, while states near the border of the is- 
land have stronger overlap with the disordered wall and 
consequently have large tunneling rates. 

Numerically, the particular realization of disorder al- 
lows for an efficient computation of quantum transport 
for remarkably long wires by employing the modular re- 
cursive Green's function method. '^^■^^ We first calculate 
the Green's functions for M = 20 rectangular modules 
with different heights. A random sequence of these mod- 
ules is connected by means of a matrix Dyson equation. 
Extremely long wires can be reached by implementing 
an "exponentiation" algorithm'*^: Instead of connecting 
the modules individually, we iteratively construct differ- 
ent generations of "supermodules" , each consisting of a 
randomly permuted sequence of M modules of the pre- 
vious generation. Repeating this process leads to the 
construction of wires with lengths that are exponentially 
increasing with the number of generations. With this 
approach we can study wires with up to 10^^ mod- 
ules, beyond which limit numerical unitarity deficien- 
cies set in. For wires with up to 10^ modules we can 
compare this supermodule technique containing pseudo- 
random sequences with truly random sequences of mod- 
ules. For configuration-averaged transport quantities the 
results are indistinguishable from each other even in the 
semiclassical limit of small de Broglie wavelength. Ad < ^, 
to be explored in the following. Accordingly, we can 
simulate quantum wires with a length of the order of 
~ 10^^ Ad, which illustrates the remarkable degree of nu- 
merical stability of the MRGM. 

The transmission tmn and reflection amplitudes Vmn 
for an electron injected from the left are evaluated by pro- 
jecting the Green's function at the Fermi energy Ep onto 
all lead modes m,n €: {1, . . . , N} in the entrance and exit 
lead, respectively. Here N = [kpW /Tr\ is the number of 
open lead modes and fcp the Fermi wave number. From 
the transmission amplitudes one obtains the dimension- 
less conductance g = Tr(ttt). The ensemble-averaged 
conductance (g) for 20 different disorder realizations and 
three neighboring values of the wavenumber fcp (Fig- 4) 
initially decreases sharply (not shown) since the contri- 
bution of modes with high m (m > 5 in the present case), 
residing primarily in the chaotic sea and transporting, on 
average, to the left, rapidly vanishes with increasing L. 

For larger lengths L the conductance shows a number 
of steps with increasing length, which can be understood 
as follows: At first, the entire regular region contributes 
to the conductance with 4 modes. With increasing L, 
fewer and fewer regular modes contribute to the conduc- 
tance, because they decay due to tunneling. This hap- 
pens first for the outermost regular mode, which has the 
largest tunneling rate jm- Finally, at very large lengths of 
the wire, even the innermost regular mode no longer con- 



tributes to the transport, so that the conductance decays 
to zero. In this limit and on this length scale, one-sided 
disorder of the quantum wire leads to localization. 

The coupling to the chaotic sea by dynamical tunnel- 
ing can also be viewed as "flooding" of the regular is- 
land by chaotic states.^'' This notion is illustrated by 
the Poincare-Husimi distributions for states injected from 
and ejected to the right (see Fig. 4, discussed below). 
The decay of the regular states into the chaotic sea or 
the reversed process, penetration of the chaotic states 
into the regular island implies that the regular states are 
quasi-discrete states embedded in a (quasi-)continuum. 
For closed wires (with periodic boundary conditions at 
X ~ and x ~ L) the discrete energy levels which form 
the quasi-continuum must feature a sufficiently high den- 
sity of states, or small level-spacing A , 

7m > A , (7) 

for effectively quenching the quasi-discrete regular state 
on the m-th torus and giving rise to chaotic eigenstates 
flooding this torus. When increasing L and thus decreas- 
ing the mean level spacing, A ~ this happens at 
characteristic lengths of the closed system where 7^ = A. 

However, Eq. 7 is not applicable to open systems. The 
step-wise behavior observed in Fig. 4 is therefore deter- 
mined by the different mode-specific tunneling rates 7™ 
themselves, rather than by a possible violation of Eq. 7. 

For the present scattering system the successive flood- 
ing of the regular island and the resulting stepwise behav- 
ior of the conductance can be visualized by considering 
averaged Poincare-Husimi distributions gained from ran- 
dom superpositions of all modes entering from the left 
and scattering to the right (insets below the curve of (g) 
in Fig. 4(a)). One clearly sees how the contribution to 
the transport from the regular region disappears with 
increasing length of the wire. Also shown are the corre- 
sponding pictures obtained for scattering from right to 
right. For small L these Poincare-Husimi functions are 
first outside the regular island which for increasing L is 
fiooded by the chaotic states. For the largest L we have 
complete flooding, i.e., no regular modes are left and the 
chaotic modes fully extend into the regular island. Thus 
the disappearance of regular states and the flooding of 
the island by chaotic states is nicely seen in the comple- 
mentarity of the Husimi pictures. 

III. LOCALIZATION LENGTHS 

We obtain quantitative results for the localization 
length ^ (in units of ^) in a wire consisting of L modules 
by analyzing the conductance g in the regime g ^ 1, 
extracting ^ from (hig) ^ —L/^. 

In addition, we can determine mode-specific localiza- 
tion lengths for the regular modes m: To this end, we 
consider the individual transmission probabilities T„i = 
J2n l^nmP for the incoming mode m as a function of L, 
see Fig. 4(b), where an average %n = exp((lnTm)) over 
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20 disorder realizations and three neighboring fc^-values 
is shown. For each of the regular modes we observe an ex- 
ponential decay subsequent to the diffusive regime. Con- 
sequently, it is possible to define the mode localization 
lengths ^„i- They can be obtained numerically by fitting 
to (In Tmm) oc —Lj^rn- The largest mode localization 
length is as the corresponding regular mode ra = \ 
couples most weakly to the chaotic region. Thus T\ de- 
termines the conductance for long wires L > such that 
^ = ^1. The regular modes with m > 1 contribute to the 
conductance for L < leading to the step- wise behav- 
ior observed in Fig. 4. 

We want to explore these localization lengths in the 
quantum-to-classical crossover. To this end we introduce 
the semiclassical parameter hcs, the ratio of Planck's con- 
stant to the area of the Poincare section, 



2pfW 



kpW 



(8) 



which coincides with the inverse number of modes. We 
study the semiclassical limit, hes ^ or equivalently 
fcp ^ oo, for two different cases: 

(i) The cyclotron radius rc is kept fixed by adjusting 
the magnetic field B = hk^/rc for increasing fcp- This 
leaves the classical dynamics invariant. In particular, the 
fraction of the regular phase space volume stays constant, 
while the absolute size Aj-cg of the regular island scales 
as Aiog oc fcp oc h~g . At the same time, Ad decreases 
and approaches the semiclassical limit. In the numerical 
computations we use the value Tc = 3W. The case m = 1 
was originally studied in our previous paper. ''^ 

(ii) The magnetic field B is kept fixed. Since the cy- 
clotron radius Vc = hk-p / B increases for increasing fcp, the 
electrons follow increasingly straight paths, thus reaching 
the disordered surface more easily. Therefore the fraction 
of the regular phase space volume decreases. As can be 
deduced from Eq. 4 in the limit ^ 0, the absolute 



area Areg still increases, Ar, 



€p (X 



-1/2 



■■eff 



but much 



slower than in the previous case. In the numerical com- 
putations we use the value B = 10.05 a. u. corresponding 
to = ?,W for h-^ = 9.6. 



A. Qualitative description 

We start with a qualitative description that gives in- 
sight in the overall dependence of the localization lengths 
for the two cases above. It is based on the observation 
that tunneling from the central mode, m = 1, of a regular 
island to the chaotic sea can be approximately described 
as being exponentially small in the ratio of island size to 
Planck's constant, ^'i^'^*^ 
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FIG. 5: (Color online) (a) Mode localization lengths for 
constant cyclotron radius, (b) Mode localization lengths 
for constant magnetic field B. The symbols connected by 
full lines display the numerical results from the full calcula- 
tion. The dashed lines show the analytical predictions from 
(a) Eq. 12 and (b) Eq. 13, while the vertical gray lines indi- 
cate the positions where the size of the regular island is large 
enough to accommodate m modes (i.e., where Arcg = mh). 



Invoking again the equivalent description of a time- 
dependent one-degree of freedom system, its tempo- 
ral decay exp(— 7it) leads to an exponential decay, 
exp(— 7ix/?;i), as a function of propagation length x. 
This gives a localization length ~ [11,12], which 
in the limit of h^s shows the following behavior. 



(i) rc fixed: exp(^co h^f^^ 



(ii) B fixed: ^ exp I cq h 
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(10) 

(11) 



with a system-dependent constant C. 



where the constant cq is different for each case and deter- 
mined below. Checking with our numerical results (full 
lines in Fig. 5) we find that, on a qualitative level, these 
estimates correctly predict an exponentially increasing 
localization length ^i. Also displayed arc the localiza- 



tion lengths for higher modes £,m, which are smaller for 
larger m. The onset of their exponential increase with 
hc« — * can be linked to the critical size of the island 
such that its area is large enough to accommodate m 
modes, A 



(ii) 



rcg 



mh. From this relation one can determine 
the corresponding values of hcs- In case (i) one finds for 
the numerically used parameters h~g w 3.5to, in case (ii) 
one finds « 0.32 + 0.59{m'^ + my/ 1.08 + m"^). These 
values arc shown in Fig. 5 by the dashed vertical lines. 

The above analysis links the exponential increase of the 
localization length to the existence of the island of regu- 
lar motion. An immediate consequence is that if such an 
island does not exist, as e.g. for two-sided disorder or for 
B = 0, an exponentially increasing localization length is 
absent.''^ For wires with two-sided disorder large localiza- 
tion lengths do, however, reappear in the quantum Hall 
regime where the cyclotron radius is much smaller than 
the wire width, Vc ^ W, and thus much smaller than 
considered throughout this paper. 



B. Quantitative description 

We now go beyond the above qualitative reasoning and 
derive analytical estimates for the localization lengths 
of the island modes. For the present realization of disor- 
der, transitions between modes only occur at the bound- 
aries between adjacent modules of differing height. At 
each boundary the wave functions to the left and to the 
right of the discrete jump have to be matched. The cor- 
responding matching conditions can be drastically sim- 
plified by considering that the lowest transverse modes in 
each module will differ only slightly from the correspond- 
ing modes in the neighboring module of different height. 
This is because the effective quadratic potential induced 
by the magnetic field, Eq. 1, plays the role of a tunnel- 
ing barrier through which only the evanescent part of the 
transverse modes may reach the upper waveguide bound- 
ary where random fluctuations of module heights occur. 
The evanescent part of the wavefunction can be repre- 
sented by a WKB-approximation for the corresponding 
tunneling integral. Within this scmiclassical description, 
we obtain an approximate analytical expression for the 
transmission coefficients and, consequently, for the local- 
ization lengths. The details of this derivation are given 
in appendix A. We find for case (i) (fixed Tc) 



(12) 



Cm ~ (am ^cff ~ cxp Co h^f^{l ~ dm hlg)^^^ 
and for case (ii) (fixed B) 

Cm ~ (am h~g^^ - brn^ CXp Co - dm h^Ji^f^'^ 

(13) 

where in each case the constants a^, bm, and dm. depend 
on m, while Cq is independent of the mode number (see 
appendix A for details). The above predictions for the 
mode localization lengths Cm arc in excellent agreement 








n 



n 



L 



FIG. 6: (Color online) Three-step process for transmission 
from regular mode m at module boundary n to the left- 
transporting chaotic modes and finally to the regular mode 
m' at the module boundary n' . 



with the numerically obtained values (Fig. 5). Moreover, 
the leading order dependence on /los is identical to the 
qualitative predictions in Eq. 10 and Eq. 11, respectively. 



IV. MODE-TO-MODE TRANSITION 
PROBABILITIES 

Additional insights into the interplay between directed 
regular motion to the right, directed chaotic motion 
to the left, and chaos-assisted tunneling can be gained 
from the mode-to-mode transition probabilities Tm'm = 
l^m'mP- These display intricate structures as a function 
of L due to the tunneling transitions between the counter- 
propagation currents the details of which wc investigate 
below. 

The diagonal transmission of a given mode into itself 
is given to lowest order by Tmm — cxp(— L/Cm)- The 
dominant contribution to the off-diagonal transmission 
probabilities Tm'm with m' ^ m, can be constructed 
based on the following three-step process (Fig. 6): (i) 
Tunneling from the right-moving regular mode m to the 
chaotic sea, (ii) propagation in the chaotic sea, which 
has an average drift to the left, and (iii) tunneling into 
the right-moving regular mode m' . This three-step pro- 
cess incorporates all the basic elements of chaos-assisted 
tunneling,^'* however here in the presence of directed reg- 
ular and chaotic transport. In addition, the signature of 
chaos-assisted tunneling is here analyzed not in terms of 
spectral properties, but in the transmission properties of 
an open systems. 



A. Qualitative analysis 

The study of mode-to-mode transmission probabilities 
as a function of the length L gives insight analogous to 
a time-dependent observation of a wave packet. Starting 
point are the transmission probabilities r,„ = J2m' '^m'm 
for each regular mode m, now plotted on a logarithmic 



scale (Fig. 7(a)), which were previously shown on a linear 
scale (Fig. 4(b)). We find that the transmission proba- 
bilities Tm which give rise to the stepwise decay of the to- 
tal transmission T also display a stepwise decay on their 
own: For short wires, the T,„ decay exponentially with 
exp(— L/^,„). Beyond L « a plateau is reached, fol- 
lowed by a slower decay exp(— L/^m-i)- This sequence 
continues until Tm finally decays with the largest local- 
ization length ^1. 

In order to explain these structures, we first consider 
the case of T4 = X^m' ^™'4 (^^^ lowest curve in Fig. 7(a)). 
Its individual contributions T^'a are shown in Fig. 7(b). 
With increasing L the dominant contribution switches 
from T44 to T34 to T24 to T14 at the localization lengths 
^4; ^3, and ^2, respectively. The appearance of the con- 
tributions T'34,T24,Ti4 can be well accounted for by the 
three-step model mentioned above (details given below). 
The subdominant plateaus occurring in Fig. 7(a) and (b) 
are due to higher order effects and will be discussed at 
the end of the section. 

The transmission of the innermost mode, Ti = 
J2m' ^m'l (shown at the top in Fig. 7(a)) and of its in- 
dividual contributions T^/i (shown in Fig. 7(c)) displays 
a qualitatively different behavior. Still, it can also be 
accounted for by the same three-step model: The trans- 
mission of mode m = 1 is always dominated by the di- 
agonal Til, as it has the largest localization length. For 
short wires the transmission probability T21 is the small- 
est of the ofF-diagonal contributions, as in step (iii) the 
mode m' = 2 couples most weakly to the chaotic sea. 
However, coupling becomes efficient with increasing wire 
length, leading to a higher saturation level than T31 and 
T41. The saturation can be understood as an equilib- 
rium of the rate of populating mode m' in step (iii) via 
the chaotic sea and the rate of depopulation due to the 
finite localization length (step (i)). The value of the 
saturation plateau will be determined from a quantitative 
analysis. 



B. Quantitative analysis 

We now demonstrate that dynamical tunneling in the 
presence of directed chaotic motion allows for a quantita- 
tive description of the transmission probabilities. While 
the opposite direction of chaotic transport in step (ii) 
was not essential for the above qualitative understand- 
ing of the transmission probabilities, it is quantitatively 
of great importance. The three steps of chaos-assisted 
tunneling (Fig. 6), can be combined to give the following 
estimate for Tm'm with m' ^ m, 

Tm'm{L) ~ X! 51 Tmm{n) — -- Tyn'm' {L - n') , (14) 

n—On'<n 

where the terms are related to the three-step process in 
the following way: 




L 

FIG. 7: (Color online) (a) Transmission probabilities Tm of 
the incoming mode m(= 1,2,3,4) vs. L showing the same 
data as in Fig. 4(b), but on a logarithmic scale, (b) Transmis- 
sion probabilities T^M- (c) Transmission probabilities T^'i- 
The dashed lines show the predictions from Eq. 14, using the 
analytical results for the localization lengths ^m. Second- 
order processes (i.e., coupling from mode m to m' and then 
from m' to n) are also taken into account. For the numeri- 
cal results, the arithmetic means (Tm) and {Tm'm) are taken. 
In regions where the distribution of transmission probabili- 
ties is log-normal, the geometric mean (i.e., arithmetic mean 
of the logarithms) could be more appropriate. For simplicity 
we use the arithmetic mean everywhere. The agreement be- 
tween analytical and numerical results is slightly better with 
arithmetic averaging. The purple ellipse indicates the region 
where the direct transition between island modes dominates 
(see text). 



(i) The first tunneling process from mode m to the 
chaotic sea can take place at any module boundary n 
along the wire consisting of L modules. The probability 
to remain in the initial mode m, before tunneling, is given 
by Trnm{n). The probability of tunneling into the chaotic 
sea is proportional to the inverse localization length . 

(ii) The chaotic sea features an average drift to the left. 
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FIG. 8: (Color online) Transmission probabilities r„'4 as in 
Fig. 7(b), but with the dashed lines showing the predictions 
from Eq. 14 with second-order processes (i.e., coupling from 
mode m to m' and then from m' to n) also taken into account. 



This is incorporated in Eq. 14 by the restriction n' < n, 
where n' is the location of the second tunneling process. 
We do not distinguish explicitly between modes of the 
chaotic sea, as they are strongly coupled and mirror the 
ergodicity of the underlying classical dynamics. 

(iii) The second tunneling process at n' from the 
chaotic sea to the regular mode m' has a probability pro- 
portional to the inverse of the localization length ^m'. 
The last factor in Eq. 14, Tm'm'{L — n'), describes the 
transmission probability within mode m' from the mod- 
ule boundary n' to the exit lead at length L. 

We now approximate the diagonal transmission prob- 
abilities Tmm by their leading-order behavior, Tmm{L) = 
exp(— L/^m). For L,^m,S,m' ^ 1 this results in 



-L 



(15) 



Eq. 15 is symmetric with respect to mode interchange, 
Tm'm{L) « T^am'iL). Three limiting cases can be de- 
duced, as 



Tm'm{L) ~ < 



■ exp {~L/^r,i) £.771' < L, ^„ 
■exp(-L/^„') £rn <.L,i,„ 



(16) 



These quantitative predictions involving chaos-assisted 
tunneling in the presence of directed chaotic motion 
are shown in Fig. 7, demonstrating excellent agreement. 
Here wc have used the analytic predictions for the local- 
ization lengths, given in the previous section. Thus, this 
agreement is achieved without any adjustable parame- 
ter. We emphasize that without explicit use of directed 
chaotic motion, e.g. without a restriction on n' in Eq. 14, 
one would get drastically different predictions not com- 
patible with the numerical results. This confirms the 
notion of chaos-assisted dynamical tunneling in the pres- 
ence of directed chaos. 



Corrections to the transmission probabilities (Eq. 15) 
can be analyzed as well. The direct tunneling between 
regular modes (i.e., without a detour to the chaotic sea) 
adds a contribution of the form 



T^'TiL) « E T7rUn)^ T^'77AL - n) , (17) 



71=0 



where the factor Q;„i'm/ (6ri$m') denotes the direct trans- 
mission probability from mode m to mode m' at a single 
module boundary. Numerical results give a factor of the 
order of am'm ~ 10, with small dependence on other 
parameters. Only for very short wires, L < 2am'7n ^ 
Cm I Cm' I Eq- 17 gives the dominant contribution to the 
transmission T^'m ~ T^J°f{L) w LQ;m'™/(CmCm'): lin- 
ear in i, which can indeed be observed in Fig. 7(c), in- 
dicated by the purple ellipse. 

Higher order tunneling processes explain the subdom- 
inant plateaus in Fig. 7(b). For example T^i{L) dis- 
plays, after the first plateau and the exponential decay 
exp(— L/Ca), further plateaus related to chaos-assisted 
tunneling from mode m = 4 to mode m" and then an- 
other chaos-assisted tunneling from mode m" to mode 
to' = 3. In general this gives for T„i'„i with to > to' > 
to" plateau values (Cm /Cm") (Cm' /Cm") in the regime 
Cm"+i <C i ^ Cm"- Fig- 8 shows the transmission prob- 
abilities T'„i'4 compared to the model predictions when 
these higher-order processes are also taken into account 
up to second order, showing excellent agreement. 



V. SUMMARY 

We have shown that in a perpendicular magnetic field 
2D nano-wires with one-sided surface disorder feature a 
regular island in phase space which leads to giant local- 
ization lengths in the limit of large Fermi momentum 
k-p , where the classical phase space structure can be fully 
resolved quantum mechanically. The coupling between 
the regular island and the chaotic sea proceeds only by 
tunneling, which is exponentially suppressed in the semi- 
classical limit. Based on this understanding, we have de- 
rived analytical results for the mode-specific localization 
length Cm in the limit of large kp , which show excellent 
agreement with the numerical data, without resorting to 
any fit parameters. 

Finally, we have presented a model describing the be- 
havior of the transmission probabilities T„m between the 
lowest modes to, n which enter and exit the wire on the 
regular island in phase space. Taking into account how 
in the interior of the wire these modes dynamically tun- 
nel to the counter-propagating chaotic sea and back to 
the island, our model shows remarkably good agreement 
with the numerical data considering its simplicity. 



APPENDIX A: ANALYTICAL ESTIMATES OF 
LOCALIZATION LENGTHS 

In this appendix we derive the analytical estimates of 
the mode localization lengths ^m- The wire consists of 
a chain of rectangular modules, with its length given in 
units of the module length I. 





1. Reduction to a single module boundary 

The transmission matrix can be constructed by con- 
necting the transmission matrices of the subsystems. 

As a building block wc first consider the connection of 
only two substructures such as, e.g. , two modules (see 
Fig. 9). The transmission matrix from left to right is 
given by 



f(2) 



.(2) 



n=0 



(Al) 



where t'-'^' , t'^' are the transmission matrices from left to 
right of the two subsystems, r'^^^ is the reflection matrix 
from the right side for system 1 and is the reflection 
matrix from the left side for system 2. 

We are interested in the transmission f^j^ of mode m 
into itself for the case that mode m is well inside the reg- 
ular island. Wc can therefore neglect all terms involving 
reflection matrices, as these would involve tunneling and 
are exponentially suppressed. This leaves 



N 

.c _/,(2),(1)n _V^(2)+(1) 
''mm ~ I- )mm — / ^ ''mn^nm 



fW +(2) 
mm mm i 



(A2) 



where N is the number of modes in the module between 
the two systems. In Eq. A2 we have neglected all but the 
mth term in that sum, as these terms involve tunneling 
to another mode and then tunneling back. 

Extending this analysis to a wire with L modules we 
get 



mm 



(A3) 



Here, tmm^^ is the transmission of the mth mode into 
itself from module i to module i + 1, with i = and 
i = L + 1 labeling the left and right lead, respectively. 
From the Onsagcr-Casimir symmetry relations'^'^^'*^ fol- 
lows that tfnm only depends on the two heights hi and 
hi^i, but not on the order in which they occur. Addi- 
tionally, for modes m well inside the regular island, only 
the exponentially suppressed tunneling tail reaches the 
upper side, so that the wave function in the module with 
larger height h can be assumed to be that of an infinitely 
high module as if unperturbed by the upper wall. This 
means that only the smaller of the two heights will influ- 
ence the transmission. Wc express the small deviation of 



FIG. 9: Connection of two scattering systems. Transmission 
proceeds by transmission through the first system {t''^^), an 
arbitrary number of reflections between the two systems (each 
of which gives a term r''^ V^-*), and then transmission through 
the second system (t'^^), leading to Eq. Al. 



tmrn^"^ from unity by the function em[^], which will be 
related to the localization lengths. 



*mm^^' = 1 - em[min(/ij, /ij+i)] 



(A4) 



Ordering the modules with increasing heights /i'^^^ ^ = 
1, M, we can rewrite Eq. A3 as a product over the M 
modules. 



M 



n(. 



M 

exp(-L^P^£™[/i('^)])^ 



(A5) 

where is the probability that at a module boundary 
the minimal height of the adjacent modules is /i^^^. 

We use again that the mth mode is exponentially sup- 
pressed at the upper boundary, from which follows that 
£m [^min] , belonging to the module with the lowest height 
/imin = h^^\ is the dominant contribution in Eq. A5. Ne- 
glecting all other contributions and using that the prob- 
ability that one of the modules at a module boundary is 
the lowest one is Pi = 2/M, t„„j is simply given by 



tr 



exp 



2L 



Ah 



m L"'min. 



(A6) 



The transmission probability Tm ~ \tmm\'^ of mode m 
decays according to Tm oc exp(— L/^m), allowing us to 
extract the mode localization length as 



M 



(A7) 



2. Transmission at the boundary of two leads 

We now calculate the reduction em[/imin] of the trans- 
mission amplitude from 1 at a module boundary with 
the lower height /imin- This can be done by calculating 
the transmission amplitude tmm for the simple system 
of two connected leads with different widths wl and wj? 
(see Fig. 10). 

We solve the Schrodinger equation in each lead and flnd 
the transmission matrix element tmm by wave function 
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FIG. 10: Two leads with widths wl and wh- 
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FIG. 11: (Color online) Effective potential V(y) and lowest 
mode in infinitely wide lead, for = 9.6 and rc = 3. In 
addition, wl = /imin for the geometry chosen in this paper is 
shown. 



matching at the boundary between the leads. Inserting 
the Hamiltonian Eq. 1 into the Schrodinger equation and 
separating the wave function as exp(ikxx)x{y) , we get 



2/ -2/0 



x(2/) = 0, 



(A8) 



which is an effective ID- Schrodinger equation for a parti- 
cle at energy i? = in the diamagnetic potential V{y) = 
EpKy - yof/rl - 1], where ya = -kx/B = -rck^/kp 
(see Fig. 11). 

The transverse mode wave functions Xm (y) , which are 
zero at the boundaries y = and y = w, satisfy the 
generalized orthogonality condition^^'**'^ 



(A9) 

where the modes Xmiy) have been normalized to carry 
unit flux, so that the scattering matrix is unitary. 

We now assume that we have an incoming (right- 
moving) wave in the mth mode in the left lead. The 
wave functions 'I'-^ in the left lead and in the right 



lead are given by 

oo 

vE'^(a;,y) = 0^+(x,y) + ^r„„0^-(x,2/), (AlO) 



vl/«(a;,y) 
where 

<l>n^{x,y) 



oo 

E 



X„ (y)cxp . 



(All) 



(A12) 



The continuity condition for the wave functions at x 
demands that 



vl/^(a:,2/) ^'^''{x,y) 



and 



dx 



ix,y) 



dx 



ix,y) 



(A13) 



(A14) 



In order to extract tmm, we multiply Eq. A13 by 
(fc^+ + 2fcFy/rc), add to it (— i) times Eq. A14, mul- 
tiply this equation by X^{y)^ ^ind integrate from y = 
to y = oo. The right hand side reduces to the general- 
ized orthogonality relation (A9) for the transverse wave 
functions in the lead and therefore simplifies to tmrm 



(A15) 



with being the generalized overlap integral between 

mode TiL in the left lead and mode n/j in the right lead, 

y 



O 



„R± 



2^Ff )x„,(2/)x„«(2/)dy, 

' C 



(A16) 

where w< is the smaller of the two lead widths wl, wr- 
For sufficiently high magnetic field B and Fermi momen- 
tum fcp , the wave functions in the mth mode are small at 
the upper boundary, so that we expect imm to be almost 
one. We define its deviation from 1 by 



= 1 - 



(A17) 



where for two leads corresponds to e,„[/ii] of the rough 
wire introduced in the previous subsection. 
Inserting Eq. A17 into the unitarity condition 



1 = 



N 

E 

n=l 



2e. 



^ ^ \tnm\ 



N 

E 



(A18) 
(A19) 



we find that the r„,„ can be at most O(^em)- The sec- 
ond term in Eq. A15 approximately corresponds to the 
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orthogonality condition (A9) for n ^ m, so that it is 
strongly suppressed. The integrals can be estimated to be 
of order 0[^/e^) by using that the difference between the 
TTith modes on the left and right side is of order 0(^e„i), 
while the left- moving modes are of order 0(1) at the up- 
per side of the wire, where X^ri'iu) differs from Xn-tiv)- 
Utilizing that the integral and r„m are both 0(^£m), the 
whole term should be of order 0(em), such that a priori it 
cannot be neglected. Numerically, we find that its mag- 
nitude does not exceed 0.22^, quite independently of the 
Fermi energy E-p and the magnetic field B. We neglect 
the second term in Eq. A15 in the following, keeping in 
mind that this will introduce an error of about 20% in 
our result for the localization lengths. 
We thus approximate 



tr. 



L.L + 



R+( 



(A20) 

which is independent of the order of the leads. Without 
loss of generality, we can choose wl < wn for the further 
calculation. Since the mth mode in the (wider) right lead 
is much less affected by the upper boundary than the mth 
mode in the (narrower) left lead, we replace the former 
by the wave function of the rnth mode in an infinitely 
wide lead, 



x^+(y) « xTiy) 



(A21) 



The transverse wave function of the mth mode in the left 
lead differs from this only very slightly, see Fig. 12, so we 
write it as 



xi^{y)-N' [xTiy) 



niy)] , (A22) 

where (Tm.{y) is negligible except near y = WL and is given 
by (T^{y) = Xm^(y) for y > wl, where Xm'^(y) = 0. 
The normalization factor N' is close to one and will be 
evaluated below. The longitudinal momentum in the left 
lead can be written as 



x.m x.ra 



(A23) 



The difference Ak^.m is very small and is numerically 
found to be 0{em)- 

Inserting (A21-A23) into Eq. A20 and extending the 
integral to y = oo (considering that the wave function 
Xni^iy) is zero for y > w^), we have 



2A:,°°+ + Afc,,,„ + 2fcF- 



N' [x^+ {y) ~ a™ (y)] x^^ (y) dy , ( A24) 



and thus 



tmm =N'[l-A + Ak^.mB] 



(A25) 










■ *\ 


















:.3^ 







FIG. 12: (Color online) Sketch of the first mode (m = 1) 
transverse wave function Xm^ (v) in ^ lead of width wl, 
Xm^iy) in an infinite lead, and the difference am{y), for 
= 9.6 and Vc — 3. The inset shows an enlargement 
of the region around y — wl- 



with 



A 



2kZt + 2fcF - hi™ (y )x™ + iy)dy, (A26) 



(A27) 



B = j [x^+iy)-cT^iy)]x^^iy)dy, 





To calculate N' , we use the flux normalization condition 
from Eq. A9, 

oo 

1 = 7V'2 y (^2k^+ + 2Afc,,™ + 2fcF^) 



[x™ + (2/)-^rn(y)]'dy, (A28) 



which we write as 

1 = 7V'2 [1 - 2A + C + 2Afc^,,„S'] , 



with 



B'= / [x^.+{y)-<J,n{y)fdy, 



oo 

C = J (^2fc-+ + 2fcF^) a^iyfdy. 



Inserting N' into Eq. A25, we obtain 

1-A + Ak^,rnB 



^1 - 2yl + 2Ak^^,nB' + C 
l-^C + Afc,,,„(B-B') , 



(A29) 

(A30) 
(A31) 

(A32) 
(A33) 
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where we used that the integrals A, B, B' ,C are much 
smaller than one. With the explicit expressions for B 
and B', the last term becomes 



Afc,,,„(B-B') = Ak,^„^ J [x^+iy)<Jm{y) - <Jm{yf] dy, 

° (A34) 
which is of higher order than integral C since cFmiy) is 
almost zero where the wave function has its maximum 
and Akx^m is already 0{em)- Dropping this term, we 
arrive at a simple expression for Em-, 



:C 



kZt + kF—] cr,n{yfdy. 



(A35) 



a. WKB approximation of cjm (y) 



To evaluate Eq. A35, we need an expression for the 
function am{y), see Fig. 12. By inserting Eq. A22 into 
Eq. A8 and imposing the boundary conditions (0) = 
Xm^(tfL) = 0, we find that am{y) has to be an eigen- 
function of the same Hamiltonian as X^rtiv)- but 
with boundary conditions (7^(0) = (Tm(oo) = 0. In- 
stead of a normalization condition, it has to satisfy 
cTmiwL) = Xm~''(^i)- Sincc the upper boundary wl is 
already deep in the classically forbidden region, we use a 
WKB approximation for our solution for (T„j (y) , 



Pm(y) 



exp ~ / P,niy')dy' y <WL 



X < 



, (A36) 



cxp - / Pm{y')dy' y>WL 



with 



we find 



[x™ + (^«l)]'x 



exp -2 / pmiy')dy' 



Pmiy)/p rn 



dy- 



''x.rn 



exp -2 / Pra{y')dy' 
y \ \ 



iy)/p 



■dy 



(A38) 

We extend the first integral to start at — oo and perform 
in the first integral the substitutions ry' = wl — y' , V — 
wl — y and r]' = y' — WL,r] = y — wl in the second. Since 
the integrals will only give a significant contribution near 
y ~ Wl, i.e., 77 = 0, we expand Pm{y) into a Taylor series 
to first order, 



Pm{wL+V) ~ Pmiwh) + p'„i[wL)ri- 



(A39) 



Defining p„ = p„i{wl), = p'^iwh) and Xm"^ = 
X^'^iwL), we find 



dTycxp (-2p„i77) X 







Tc J I- WrnlPrn 
WL+V\ {-p'rnV^) 



,00+ 



1 + Wm/Pr, 



(A40) 



Expanding the term in square brackets in powers of i] 
gives 



kpWL\ 
Tc J 



0{if). 



(A41) 



Dropping the quadratic term, the evaluation of the inte- 
gral leads to 



(x: 



oo+\ 



00+ 



k-pWL 



Pm \ 'cj 

with Xm'^ and pm both evaluated at y — w^. 



(A42) 



Pm {y) = V2[Kn(y)-S] = fcp ^/VM/Ep, {Ail) 

where we have used that the effective ID Schrodinger 
equation has eigenvalue E = Q and the effective potential 
Vra depends on the longitudinal wave number fc^.m (cf. 
Eq. A8). Inserting the above expression into Eq. A35, 



b. WKB approximation of modes Xm^iv) 

The next step in our calculation of £„, and, ultimately, 
of the mode localization length is to find an expres- 
sion for the value Xm~'"(^i) of the transverse wave func- 
tion at the upper boundary and for the longitudinal mo- 
mentum eigenvalue k"^^. 
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To do this, we rewrite the effective transverse Hamil- 
tonian Ry , following from Eq. A8, 



small error. Therefore, we insert Eq. A47 into the flux 
normalization condition Eq. A9 and obtain 



with the potential 



V{y) 



y-Vz 



(A43) 



(A44) 



with the classical turning point j/z = yo + Vc, where the 
potential value V{yz) = equals the energy E = 0. We 
linearize the potential near the classical turning point 
(where the WKB solution diverges), which makes it pos- 
sible to solve the efi^ective Schrodinger equation analyti- 
cally. The complete wave function is constructed by con- 
necting the solution of the linearized potential near the 
classical turning point to the WKB solution in the clas- 
sically forbidden region. 

The linearized Schrodinger equation is solved by the 
Airy function Ai(z) with z = (2kp/rc)^^^{y — y^). 
The boundary condition at the lower wall demands 
that x{y = 0) is equal to zero, so that z{y = 0) 
must be a zero of Ai{z). For the solution for the 
mth mode, we choose the mth zero at z = w 
(-2.338,-4.088,-5.521,-6.787,...), so that the wave 
function has m — 1 nodes. From this follows that 



2fcp 



1/3 



and by using yo 



ckx/k-p we find 



1 + ^ ('^F^-) 



(A45) 



(A46) 



For large k-p the longitudinal momentum of the mth mode 
k'^^ is only marginally smaller than the Fermi momen- 
tum k-p (remember that z„i is negative). The maximum 
of the mth transverse mode is between z = Zm and z = 0, 
i.e., between y = and y ~ yz- Since approaches y — 
for large fcp, the wave function stays closer and closer to 
the lower wall with increasing kp. 

The transverse wave function near the classical turning 
point can now be written as 



1/3 



(A47) 



We will construct the full solution for the wave func- 
tion by using the Airy function near the classical turning 
point and the WKB solution (which takes the quadratic 
potential into account) in the classically forbidden re- 
gion. Before proceeding, we determine the prefactor Cm- 
Since the WKB solution is only used in describing the ex- 
ponential tail for y ^ yz, calculating Cm with the wave 
function of the linearized potential will only introduce a 



oo 

1 = j (2k^+ + 2kp^^ [x5^+(y)] ' dy, (A48) 



which results in 



Cn 



{Akpr, 



Nl/3 



21/3 {kpr,f'\ 



-1/2 



(A49) 

oo oo 

where /i„i ~ J Ai^(z)dz and i^m = J z Ai^(z)dz. In the 

Zm Zjn 

limit of large kp, this simplifies to 

-1/2 



c„ 



(A50) 



Since we need to evaluate the transverse wave function 
Xm^ at 2/ = , which is deep in the classically forbidden 
region, we proceed by connecting the Airy function (valid 
near the classical turning point) to the WKB solution 
(valid in the classically forbidden region). We write the 
WKB solution as 



xTiy) 



V 

:cxp I - J Pm{y')dy' 



y > yz, 



(A51) 

and from a short calculation we obtain that the two con- 
stants Cm and Dm are related by 



2kl 



2 \ 1/6 



Crr 



20F' 



(A52) 



For evaluating the integral in Eq. A51, we insert the 
explicit form of the potential, use y = wl as the upper 
limit of integration, and rewrite Eq. A37 as 



Pmiy)^kpd2——+l—— 



(A53) 



leading to 



z 

Pmiy'W ^kpr, j V2z' + z'Mz', 



(A54) 



with z = (wl — yz)l'^c ■ We evaluate the above integral 
by expanding the integrand in powers of z'. 



1 PmWy' - V^fcpre 1 (^Z'^/' + ^ + 0(z^/2) j dz' 



2^2, 



kprj^l'^ ' ^ + 20^ 



(A55) 
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Inserting this into Eq. A51, and using that <C wl 
in the hmit of large k-p, we obtain 



Dr 



\/PmXwL) 



exp 



2V2 



, 3/2 / X 3/2 

Wl 



1 + 1^ 

20 Tr 



(A56) 



ncglceting higher order terms. 

We now insert this into Eq. A42 and obtain 



Em exp 



4V2, (wL 



20 



3/2 



Ml. 

WL 



3/2' 



with the prefactor 



(A57) 



(A58) 



21/3471^™ [wL - V^f/rl - 1 ' 
where we used Eq. A52, Eq. A50, and Eq. A37. 



3. Semiclassical limit with constant cyclotron 
radius 

To obtain our final result for the mode localization 
length S^rn in the semiclassical limit of large k-p , we insert 
Eq. A57 with wl = /imin into the expression Eq. A7 for 
^rn- We keep the cyclotron radius Tc = ck-p/B constant, 
such that the classical dynamics is independent of kp. 
We expand the prefactor M/{AEm) in powers of kp for 
/cf ^ 00, keeping the first two terms in the expansion 
since they are of similar magnitude for the parameter 
values used, and finally obtain 

(A59) 

where h^g = (kpW/'!r)~^ and the dimcnsionlcss parame- 
ters flm, 6m, Co, and dm arc given by 



2 + A/C \ 



ci/3 vi+A/cy 



Co 



4^/27r 



— - 

3 C^/^ V ^ 20 C 



c 



1/3 



21/3,^2/3 A 



(A60a) 
(A60b) 
(A60c) 
(A60d) 



where we have introduced the dimensionless parameters 
A = /imm/VK = l-{5/2W) and C = rjW. Note that the 
factor Co determining the asymptotic exponential behav- 
ior does not depend on the mode number m. For m = 1 
and in the limit A/C <C 1 this corresponds to the result 
we have previously reported. 



4. Semiclassical limit with constant magnetic field 

Instead of keeping the cyclotron radius fixed, we al- 
ternatively set the magnetic field i? to a fixed value in- 
dependent of kp and again perform the limit kp 00. 
We repeat the above procedure, introducing the magnetic 
length A through B/c= 1/A^ and the dimensionless pa- 
rameter A = \/W . As above, we use A = /imin/W^ and 
obtain 



6, 



7-1/3 _ , 



,)cxp(co (l-d™/i^/')3/2), 

(A61) 

with the dimensionless parameters am, bm, cq, dm now 
given by 



am - {2TTf^Hl^l„ 
bm = ~2TTZmMfJ,ri 

4^/2^ A3/2 



A2/3' 



Co 



A ' 

^2/3 



(27r)i/3 A 



(A62a) 
(A62b) 

(A62c) 
(A62d) 
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